Spin-dependent transport in a Luttinger liquid 
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We develop a detailed theory for spin transport in a one-dimensional quantum wire described 
by Luttinger liquid theory. A hydrodynamic description for the quantum wire is supplemented by 
boundary conditions taking into account the exchange coupling between the magnetization of ferro- 
magnetic reservoirs and the boundary magnetization in the wire. Spin-charge separation is shown 
to imply drastic and qualitative consequences for spin-dependent transport. In particular, the spin 
accumulation effect is quenched except for fine-tuned parameter regimes. We propose several feasi- 
ble setups involving an external magnetic field to detect this phenomenon in transport experiments 
on single-wall carbon nanotubes. In addition, electron-electron backscattering processes, which do 
not have an important effect on thermodynamic properties or charge transport, are shown to modify 
spin-dependent transport through long quantum wires in a crucial way. 

PACS: 71.10.Pm, 72.10.-d, 75.70.Pa 



I. INTRODUCTION 

Spin-polarized transport represents a new branch of 
mesoscopic physics, in which both the charge and spin 
of the electron are actively manipulated. In differ- 
ent setups, the spin may be employed as information 
storage or transport element, where the advantage over 
charge transport stems from the very long spin lifetimes 
in many materials, and the smallness of the dissipated 
power. These advantages have already resulted in many 
technological applications, and among the most popu- 
lar future perpectiyes of spintronics is the field of quan- 
tum computation.!! Spin-dependent transport also offers 
novel insights into fundamental physics. In this paper 
we shall address in detail how spin transport proceeds in 
strongly interacting non-Fermi liquid metals, taking the 
behavior of one-dimensional (ID) metals as a paradigm in 
which electron-electron interactions lead to a breakdown 
of Fermi liquid theory. The ID non- Fermi liquid behavior 
is often described by Luttinger liquid (LL) theory.0 

The primary motivation far this study comes from 
recent transport experimentsu on carbon nanotubes ,u 
which have demonstrated the breakdown of Fermi liquid 
theory in these nearly ideal ID quantum wires (QWs). 
In fact, when studying charge transport in single- wall 
nanotubes (SWNTs), the observed power-law behaviors 
in the tunneling density ofj-states are consistent with 
their theoretical description^ in terms of a LL. The 
LL describes metals in the ID limit where only one or 
very few bands intersect the Fermi energy. This non- 
Fermi liquid exhibits fractionalization of electrons into 
novel quasiparticles, comprising a diverse set carrying 
spin separately from .charge, and charge in fractions of 
the electron charge e.Eil Furthermore, the LL is the sim- 
plest model showing the remarkable phenomenon of spin- 
charge separation which has been postulated, by many 
to under ly the cuprate superconductors .tHnlj In a LL, 



spin and charge degrees of freedom are completely decou- 
pled, and moreover characterized by different velocities. 
As Landau quasiparticles are unstable, an electron will 
spontaneously decay into charge and spin density wave 
packets which then propagate with different velocities. 
Thereby a spatial separation of spin and charge of the 
electron results. Unfortunately, this hallmark behavior of 
strongly correlated ID fermions remains to be observed 
experimentally (at least in an unambiguously accepted 
way). Elaborating on our short paper ,E3 we propose sev- 
eral feasible setups that would allow to unambiguously 
detect spin-charge separation via spin transport exper- 
iments on an individual SWNT. For related but rather 
different proposals to detect spin-charge separation in a 
LL via spin transport, see also Refs. [jj] and [l8| 

As will be described at length below, this can be 
achieved by attaching ferromagnetic leads to the QW, 
and possibly an additional magnetic field, see Fig. [|. 
For simplicity, identical contact parameters are assumed 
below, with straightforward generalization possible. By 
measuring the variations of the current-voltage (/ — V) 
characteristics with either the angle 9 between the fer- 
romagnetic magnetizations in the leads, or the magnetic 
field B, one can indeed directly probe spin-charge sep- 
aration. A related spin-transport experiment has been 
carried out recently for multi-wall nanotubes, where the 
angle 9 was fixed to either zero or 9 = 7r.t£l The exper- 
iment proposed here for SWNTs should either allow for 
arbitrary 9, or employ an additional magnetic field. We 
note that spin transport in such a setup is well under- 
stood for Fermi liquids. In particular, the I — V charac- 
teristics (including a magnetic field) have .been computed 
recently using a semiclassical description. E3 

In our theory, we assume that tunneling across the two 
contacts proceeds incoherently, i.e., the length L of the 
QW must be longer than either the thermal length scale 
hv/ksT or the scale fiv/eV set by the applied voltage 
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V. We then consider in general systems composed of ID 
interacting quantum wires and bulk ferromagnets. With 
the exception of Sec. VI, where the additional flavor de- 



gree of freedom present in SWNTs will be addressed, we 
focus on single-mode QWs. 

Another interesting aspect of spin transport concerns 
the role of electron-electron backscattering interactions. 
In a spin-i QW, these interactions are (marginally) ir- 
relevant under the renormalization group (RG) flow, and 
therefore only cause a renormalization of interaction pa- 
rameters in the low-energy LL theoryE However, this 
essentially thermodynamic (equilibrium) argument must 
be re-examined when dealing with spin transport. In 
fact, such interactions, despite being irrelevant, can re- 
sult in nonlinear and sometimes dramatic effects, e.g., a 
non-sinusoidal oscillatory I — V characteristics. 

The structure of this paper is as follows. In Sec. [Sj, we 
introduce the basic model and outline the computation of 
the non-equilibrium spin current. In Sec. Ill, the physics 
arising at a contact between a ferromagnetic reservoir 
and a Luttinger liquid is addressed at length. Two pro- 
cesses are shown to be of importance, namely electron 
tunneling and boundary exchange. Exchange leads to 
novel conformally-invariant boundary conditions which 
are derived here. In Sec. IV, a hydrodynamic description 
of spin transport in the ID QW is developed. In Sec. 0, 
we derive the I—V characteristics in a magnetic field for 
the simplest spin transport device, see Fig. [I], first for a 
short-to-intermediate length of the QW. Under the latter 
condition, backscattering can be neglected. The effec ts of 
bac kscat tering are then addressed in detail in Sees. VC 
and VD, where we focus on zero magnetic field for clar- 
ity. Finally, several extensions and possible conc erns are 
addressed in Sec. VI. We conclude in Sec. VII by dis- 



cussing an analogy to ballistic superconductor-normal- 
superconductor (SNS) junctions, summarizing some open 
questions and providing an outlook. Details of our cal- 
culations in Sec. ^ can be found in three appendices. In 
intermediate steps of the calculations, we put e = % = 1, 
but restore units in experimentally relevant results. 



II. MODEL AND FORMULATION 



The low-energy description of a single-mode QW is re- 
markably universal ,0 and a sufficient Hamiltonian for our 
purposes is 



^-i^ VT z d x ip + u(^ip) 2 ^ , (2.1) 



Hqw = I dx 



where ip = ip aa is a four-component spinor. Here a = 
R/L indexes the chirality that differentiates right- and 
left-moving modes, and a =f / J, indexes the spin. Wc 
suppress the indices whenever possible, employing Pauli 
matrices f and a acting in the chirality and spin spaces, 
respectively. Furthermore, v denotes the Fermi velocity, 
or, more generally, the spin velocity of the interacting 



theory. If the QW is isolated, the zero-current condition 
at the endpoints requires to impose the boundary con- 
ditions ipn(0) — "0l(O) and iPr(L) — iPl{L). Only the 
forward-scattering interaction u is kept in Eq. ( j2.l| ). Al- 
ternatively, we will use the exponent a > for tunneling 
into the end of the LL, e.g., at x = 0, as a measure of 



the interaction strength. Eq. (2.1) is not completely gen- 
eral. It contains two parameters, v and u, while a general 
single-mode LL has three parameters: a charge velocity 
v c , a spin velocity v, and a dimensionless "Luttinger pa- 
rameter" ,often denoted K p or g, where a = {K'Z 1 — 1) /2. 



In Eq. (2.1), we have assumed full (Galilean) transla- 
tional invariance, leading to K p = v/v c . This relation 
is expected to be well-satisfied in many experimentally 
relevant QWs, and moreover the manipulations to follow 
relax this condition and thus can be applied even when 
K p ^ v/v c . 



In some circumstances, Eq. (2.1) should be supple- 
mented by the electron-electron backscattering interac- 
tion, 



with the chiral spin currents 



bv / dx Jl ■ Jr 



(2.2) 
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Jr/l(x) = I : *l) R/L {x)3ip R/L {x) : , (2.3) 

where the colons denote normal ordering. These chi- 
ral spin currents obey Kac-Moody commutation relations 

(H,v = x,y,z), 

+ ie^ x Jt /R {x)5{x - x') , 

where the + (— ) sign is associated with the L (R) current. 
We not e th at in a spin- 1/2 QW, the backscattering inter- 
action ( |2.2| ) is marginally irrelevant in the RG sense, and 
hence can be neglected at low energies in many equilib- 
rium properties. In a SWNT, the generalization of ifbs 
causes exponentially small gaps that can be neglected 
at not too low energies. Furthermore, the dimensionless 
backscattering coupling constant b is generally small and 
scales as 1 /R with the tube radius R. Importantly, as will 
be discussed below in detail, a precession effect encoded 
in Eq. (^^) is crucial for understanding spin transport in 
long QWs. 

For energies well below the electronic bandwidth D, a 
ferromagnetic (FM) lead can be describecUising an effec- 
tively non-interacting Stoner-like picturecil with a con- 
stant density of states. It is then sufficient to employ a 
non-interacting ID model, e.g., for the left lead (x < 0), 



Hfm = 



E 

s=±l 



dxf (-iv s u s T z d x ) f , 



(2.5) 



where / is again a four-component spinor. Comparing 
to Eq. (2.1), different spin quantization axes have been 
used, and therefore the projection operator 
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u s = (l±m-<x)/2 



(2.6) 



M = J R + J L 



(2.11) 



projecting the spin quantization axis of the QW onto 
the magnetization rh is needed. This description of the 
semi-infinite lead must be supplemented by an appropri- 
ate boundary condition, Jr(0) — /l(0). In Eq. (2.5), 
the two Fermi velocities v± parameterize the different 
densities of states, p a = l/(2nv s ), for the majority and 
minority carriers. Following Ref. |2^, we choose a suit- 
able rescaling of the / operators to set v + = V- = 1, 
thereby incorporating the difference in the density of 
states into a redefinition of the hopping matrix elements, 
t s —> t s / p s , employed in the tunneling Hamiltonian, see 
Eq. (2.7) below. Formally this is done by choosing eigen- 
states f s of rh ■ & with eigenvalue s = ±, and then rescal- 
ing f s (x) —> fg(x/v 8 ), the spatial rescaling being al- 
lowed because the different spin polarizations are non- 
interacting and tunneling acts only at x = 0. 

The LLs and FMs in question will be considered cou- 
pled by low- conductance contacts with identical proper- 
ties. Processes in which electrons are transferred across 
such a contact can be described by the tunneling Hamil- 
tonian -fftun- Provided this contact occurs at one of the 
ends of the LL, say at x = 0, this has the form 



H tun = F^W^ + VWF , 



(2.7) 



where F = /(CP) and if? = ip(0 + ) are Fermion annihila- 
tion operators at the ends of the FM and QW, respec- 
tively. The 2x2 tunneling matrix W reads with Eq. (2.6) 



W 



E 

s=±l 



t s U s 



(2.8) 



Using the spin-dependent hopping matrix elements t s , 
we may define spin-dependent conductances G s = 
(e 2 /h)\t s \ 2 , or, alternatively, the contact parameters 



G = G T +G X , P = (G T -G X )/G 



(2.9) 



Here G is the total conductance associated with the con- 
tact. In a slight abuse of terminology, we call P the po- 
larization. The polarization satisfies < P < 1, and in 
fact represents the asymmetry between the local tunnel- 
ing density of states of the majority and minority carriers 
of the ferromagnet. Hence P is not a bulk property, and 
depends upon the detailed nature of the FM-QW con- 
tact. In the experiment of Ref. |l9|, application of our 
theoretical results, in particular Eq. ( |5.9|) - which differs 
slightly from the theory used in Ref. |19| - gives a value 
of P = 0.3 for a multi-wall NT to FM (cobalt) contact 
at T = A.2K. 



From Eq. (2/7), one may deduce the spin current, de- 
fined by 



where the magnetization density in the QW is 



(2.10) 



Using the continuity equation in the absence of backscat- 
tering interactions, 6 = 0, the steady-state current is thus 



J = v(J R - J L ) . 
Then the tunneling current is 

J tun = 1 {F*W&U - &ffW*F) 



(2.12) 



(2.13) 



Of course, a formula similar to Eq. ( ]2 . 13| ) obtains for the 
charge current across the contact, 



I = i (F^W^ - yiW^F) 



(2.14) 



In addition to the tunnel coupling in Eq. (2.7), the FM 



magnetization and the LL boundary magnetization can 
be coupled by a pure exchange term, 



-Km ■ VP 7 -* 
2 



(2.15) 



Even if in a microscopic formulation no bare exchange 
coupling K is present, it will be generated in the low- 
energy effective Hamiltonian since tunneling causes vir- 
tual processes corresponding to exchange, see below and 
Ref. ||. 

To study transport, we must formulate the non- 
equilibrium dynamics of the system. This formulation 
is more subtle than in conventional charge transport, 
due to the complications arising from the non-commuting 
nature of spin. We therefore proceed carefully along 
the lines of a Keldysh approach. For concreteness, we 
specialize for the moment to a semi-infinite LL con- 
tacted at x — 0. Consider an initial system composed 
of two decoupled pieces, described by the Hamiltonian 
Hq = Hfm + Hqw- The ferromagnetic part, governed 
by ifpM, is polarized along direction rh, and located at 
x < 0. The wire, located at x > 0, is governed by Hqw, 
which is assumed SU{2) invariant. The latter condi- 
tion guarantees the existence of a continuity equation for 
the spin density and current. Similarly, charge conser- 
vation implies a continuity equation for charge density 
and current. At t — — oo, we assume that each half is 
at quasi-equilibrium at its own chemical potential, /xfm 
and (J.QW- Similarly, we assume that the LL supports a 
quasi-equilibrium magnetization, which can be described 
by a grand canonical distribution with a "spin chemical 
potential" h. We stress that neither /xqw n or h are physi- 
cal potentials, such as electrostatic or Zeeman fields, but 
rather characterize the initial non-equilibrium distribu- 
tion. Then we adiabatically turn on the contact pertur- 
bation H' = H tun + H cx , 



H{t) = H + e 5t H' , 



(2.16) 



where S — > + is an infinitesimal inverse time-scale con- 
trolling the slow turning on of the contact interaction. 
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The above formulation is rather rigorous, but has the 
limitation of being formulated to treat both the tunnel- 
ing and boundary exchange terms as perturbations. In 
practice, this is appropriate in the low- conductance limit, 
since the boundary exchange is generally determined by 
virtual tunneling processes, and is hence small. However, 
it is useful theoretically to contemplate a situation in 
which the tunneling is small but the boundary exchange 
is not. In such a case, it appears natural to include H cx 
into Hq rather than in H'. This raises difficult concep- 
tual issues, since H ex does not commute with the total 
spin of the LL, and hence renders its magnetization un- 
certain in directions perpendicular to to. On physical 
grounds, however, we expect that, because H ex transfers 
spin into/out of the LL only at the boundary, sense can 
be made of the bulk spin chemical potential in the limit 
of a semi-infinite LL. This is fairly clear from the follow- 
ing thought experiment: imagine preparing the LL with 
zero tunneling and zero boundary exchange in a state 
with a non-zero magnetization density not parallel to to. 
Then, if H cx is turned on at some time, its effect will be to 
scatter left-moving electrons into right-moving ones upon 
their reaching x — 0, changing their spin orientation in 
the process. If the LL is semi-infinite, however, it would 
take an infinite amount of time to modify the mean mag- 
netization of the LL in this way. Instead, one expects 
a steady-state to be established, generally with a time- 
independent spin current. Moreover, for a finite but long 
LL, so long as there are some inelastic processes deep 
in the LL that can equilibrate the returning electrons, 
one expects that an equilibrium state will be established 
which has a mean magnetization very close to that before 
turning on the boundary exchange. In fact, a scattering 
approach of this type can be directly implemented using 
bosonization methods, and will be discussed analytically 
in Sec. 



IIIB 



Alternatively, H ex can be incorporated di- 
rectly into Hq, but in this case, care must be taken to 
ensure that h is coupled only to the magnetization out- 
side a neighborhood of the boundary. For the moment, 
we simply use the above discussion as motivation to in- 
corporate Hex into Hq. 

We are interested in the properties of the system at 
time t = 0, when a steady-state transport of charge and 
spin between the two systems has been achieved. Then 
we can formally calculate the expectation value of any 
operator at this time: 



(O(0)) = ^Tr 



-P(H - 



A >Tc 



xTexp 



exp 



dtH(t) 



dtH(t) ) O 



(2.17) 



(2.18) 



where T denotes time ordering and 

A = hfmNfm + MQw^qw + h ■ Sq^j . 

It will be most convenient to choose the zero of energy 
in H such that /ifm = Mqw = h = in equilibrium, i.e., 
for zero applied voltage. 



Expanding out the time-ordered exponential, one ob- 
tains to lowest order 



(O) =~i dt 



iTr (e-P"°[0,H'(t)} 



(2.19) 



where Hq = Hq - A and H'{t) = e lHot H' e~ iHat . Owing 
to the factor of A in the Boltzmann weight, not present 
in the time evolution of H'(t), this is a non-equilibrium 
expectation value. We can cast it into a more equilibrium 
form by writing Hq — Hq + A. Then 



H'(t) = e l " ot (e lAt H'e- lAt ) 



-Hint 



(2.20) 



Equation (2.2C) has the form of the time-evolution 



of the operator in the brackets evolved by the "non- 
equilibrium" Hamiltonian H 0l which i s the same as is 
used in the Boltzmann weight in Eq. ( [2.19 ). Equation 
( |2.19 ) can be rewritten in the slightly more suggestive 
form, 



(O) 



dte At ([0,H' A (t)]) A 



Ho 



(2.21) 



where H' A indicates the modified operator 



H' A = e lAt H'e~ lAt 



(2.22) 



which evolves according_to the fictitious equilibrium time 
evolution dictated by Hq. The subscript Hq on the ex- 
pectation value indicates that it is a standard equilibrium 
average with respect to the Hamiltonian Hq, in which the 
argument of the Fermion fields indicates standard Heisen- 
berg picture time dependence using Hq, 
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(Oi'-O n ) 6a = -Tr 



-PHo 



Ox ■ ■ ■ O t 



z 

0(t) = e i6at Oe- i&ot . 

Thereby we can express an intrinsically non-equilibrium 
property of the system with Hamiltonian Hq in terms of a 
fictitious equilibrium average with respect to the shifted 
Hamiltonian Hq. 



III. CONTACTS 



In this section, we analyze in detail the physics of a 
single contact between a FM lead and a semi-infinite LL 
(taken at x = 0). We expect that the results apply to 
finite-length LLs longer than the thermal length v/ksT 
beyond which transport is incoherent. 
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A. End contacts: Boundary operators 

We first study the properties of the contact in equilib- 
rium from the RG point of view. It is helpful to view 
both -f/tun and H cx as perturbations to a decoupled fixed 
point described by Hq- Standard arguments give the scal- 
ing dimension of both t s and K , 



l + a/2, A K = l 



(3.1) 



The scaling dimension A^- is not renormalized due to 
spin-charge separation in the QW. A simple calculation 
then gives the RG scaling equations 



d e \t s \ 2 (£) = -a\t s \< 



d l K(e)^c(\t ] f-\t l \ 2 



(3.2) 



where £ — \n(D/E) is the standard RG flow parameter, 
and c denotes a non- universal constant. Note that K is 
renormalized by the hopping matrix elements, but the 
converse does not occur. The most important property 
of Eq. ([T^) is that while the tunneling is irrelevant for 
a > 0, the exchange coupling K between the LL and FM 
is exactly marginal. 

Following the RG flow from the ultraviolet cutoff D 
down to energy E max(fcgT, eV) <C D, we find 
\t s \ 2 (E) = \t s \ 2 (E/D) a , and 

K{E) =K + ar l cGP[l - (E/D) a ] > \t s \ 2 (E) . (3.3) 

Therefore the effective exchange coupling K(E) does not 
pick up the (E/D) a suppression factor, and generally 
is much larger than the effective hopping t s (E), regard- 
less of the microscopic "bare" values of these couplings. 
Parenthetically, we note that the non-interacting limit, 
a — > 0, is not correctly handled by the simple RG equa- 



tions (3.2), as this limit involves additional marginal op- 
erators. In fact, for a — > 0, not logarithmic dependencies 



[as predicted by Eq. (3.3)] but instead a principal parts 
prescription emerges. 

Since the exchange coupling is exactly marginal and 
the tunneling irrelevant, solving the problem for zero 
tunneling but non-zero exchange gives the "boundary 
fixed point" solution. From the viewpoint of low-energy 
physics, the only effect of the boundary exchange cou- 
pling is then to induce a modified boundary condition at 
the contact. This modified conformally invariant bound-, 
ary condition comprises a novel boundary fixed pointed 
describing the semi-infinite LL close to a ferromagnet. 



B. Zero Tunneling Boundary Fixed Point 

To gain maximum insight into the physics, we solve 
the equilibrium problem with zero tunneling exactly in 
a number of equivalent ways. The most familiar method 
is Abelian bosonization.Q Choosing the quantization axis 



for the spinor basis along the m axis, the electron field 
can be written in terms of boson fields, 



where (ip p ,0 p ) and (ip a ,6„) are charge and spin bosons, 
respectively, that satisfy the algebra 

with the Heaviside step function <d(x). The short- 
distance cutoff ao describes a non-universal scale factor 
relating the microscopic Fermion field to the continuum 
bosonized vertex operators and is related to the band- 
width, clq ~ v/D. 



Neglecting the bulk backscattering interaction (2.2) for 
the moment, the LL Hamiltonian splits into spin and 
charge components, H = H p + H a . We require only the 
spin component, 



Ha= 2^J dx ^ dx9 ^ 2 



K 
ttV2 



d x d a (0) . 

(3.5) 



The last term is the boundary exchange term. It can be 
transformed away by the canonical transformation 



r (x) ^ 6 a (x) + -^=Q(x - 0+) 



(3.6) 



which simultaneously encapsulates several physical ef- 
fects. First, since 

S-m~d x 6 a /{-KV2) , 

a local "proximity effect" magnetization is induced in 
the neighborhood of the contact. As this induced mag- 
netization decays on the microscopic scale of the Fermi 
wavelength, the magnetization appears within bosoniza- 
tion as a delta function. Second, the transverse left- and 
right-moving spin currents also depend on 9 a , 



J 



Thus the shift in Eq. (3.6) leads to modifications of the 
transverse spin current, 

J!(x>0+)^e ± " 5 / 2 j!(x>0+), (3.7) 
J±(x>0+)^e TJ,?/2 J±(x>0 + ) , (3.8) 

where a dimensionless measure of the boundary exchange 
coupling is provided by the exchange angle, 



<d = 2K/v 



(3.9) 



This transformation can thus be interpreted physically 
as a phase shift. Where before the transformation spin 
conservation required Jr(0 + ) = Jl(0 + ), after the change 
of variables, we have 
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</|(o + ) 



(3.10) 



Unlike the purely local magnetization parallel to rh, this 
phase shift can have measurable consequences far from 
the contact. 

The phase shift can also be understood directly in 
terms of electrons, which is useful far making contact 
with earlier non-interacting theories.E3 It is simplest to 
combine the right- and left-moving Fermions of the semi- 
infinite LL into a single chiral right-moving Fermion for 
each spin species on the full infinite line. In particular, 
we let 



tpLa{-x) 



(x < 0) 
(x > 0) 



(3.11) 



which ensures continuity of ip' at the origin due to the 
boundary condition tpR(0) = ^l(O). At the bound- 
ary where the original right- and left-movers have been 
"merged" together, we have Jr(0) = Jl(0) = J'(0), so 
that the exchange Hamiltonian becomes 



H cx = -2Km ■ J'(x = 0) 



(3.12) 



with J' = (tp'y (Tip' /2. We may then write down the 
Dirac equation for ip' as 



(<9 f + vd x )tp'(x) — iKS(x)rh ■ atp'(x) 



(3.13) 



For low-energy station ary st ates, the time derivative may 
be neglected, and Eq. ( 3.13j ) then yields 



ijj'{x = 0+) = cxp U-m • aj ijj'{x = 0") . (3.14) 

Thus the boundary exchange simply induces different 
phase shifts for (left-moving) electrons incident upon the 
contact from the LL and reflected back into the LL (as 
right-movers) , dependent upon their polarization relative 
to rh. 

It is technically most convenient to work directly with 
the spin currents. This has the advantage of keeping the 
spin quantization axis arbitrary at all stages. Using the 
same "merged" operato rs d efined abo ve, th e Kac-Moody 
commutation relations (2^4) and Eq. ( 3.12 ) result in the 
equation of motion for the merged chiral spin current, 



(d t + vd x )J' = 2KS(x)m x J' 



(3.15) 



where bulk backscattering is again neglected. In a steady 
state, d t J' = 0, and Eq. (|j5|) can be formally solved to 
obtain 



MO + )=K(0)J L {O + ) 



(3.16) 



Here the phase shift is encoded in the one-parameter 
50(3) matrix, 

H{&) = exp(tfF) , IV = "AfAp , (3-17) 

A 



which describes rotatio n by an an gle d around the rota- 
tion axis rh. Equations (3.16) and ( p. 17 ) provide the most 
general formulation of the effects of boundary exchange. 
Using them, we take i? to define the dimensionless "ex- 
change coupling constant" of the low-energy theory. It 
describes the angle that an incident spin in the LL pre- 
cesses around the FM magnetization direction rh due to 
the exchange interaction. In principle, since the bound- 
ary exchange operator is exactly marginal, d need not be 
small, but for the case of low-conductance contacts, one 
has •& •< 1, see Ref. ^2[ The boundary condition ( 3.16j ) 
describes the zero tunneling boundary fixed point in the 
presence of exchange and is crucial to the subsequent de- 
velopment of our theory. 



C. General formulation 

We now include the effect of tunneling on top of the 
boundary exchange above. To do this is slightly subtle, 
owing to an (unphysical) short-distance singularity in- 
herent to the linearized spectrum of the Luttinger model. 
To resolve this difficulty, we are required to choose some 
short-distance regularization for the microscopic physics 
of the contact. The form of the resulting macroscopic 
equations is independent of this choice, although the 
quantitative values of certain 0(1) coefficients can be 
cut-off dependent. A convenient method is to employ the 
combined infinite chiral Fermion description introduced 
above, and then assume that the tunneling occurs on the 
right-moving branch slightly after the exchange coupling 
acts, i.e., within some distance of the order of a,Q. When 
an electron tunnels into the LL from the FM, its spin and 
charge are propagated to the right and do not themselves 
interact with the exchange torque, so that 

Jl = Jl{0 + ), Jr = 4(0+) + - J tun ■ (3.18) 

V 



The additional tunneling spin current, Jtun, can now be 
calculated using the time-dependent perturbation theory 
treatment described in Sec. |HJ. 

In particular, we consider O = Jt un wi th Eq. ( 2.13 ), 
and H' = H tnn . For this case, Eq. (2.22) yields H' A = 
A(t) + A*(t), where 



A(t) = F^WU{t)m , 
with the unitary matrix 



U(t) = cxp 



V 



(3.19) 



where V — /Uqw — A'fm- Note that U(t) is simply a 
matrix, and hence does not represent an operator in 
the Hilbert space. It comprises the only explicit time- 
dependence in the integrand in Eq. ( 2.2l| ), and can be 
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remo ved o utside the trace. Applying the above results to 
Eq. ( 2.21 ), we find (repeated indices are summed) 



(J tun ) = Re / dte st (Wa) a0 (U\t)W*) 



7 A 



Fl(0)VJO)M(t),FJt) 



H 



(3.20) 



A formu la si milar to Eq. ( 3.20 ) obtains for the charge 
current ( [2.14 ), 



(/) = 2Rc / dte 5t W aP (rf{t)W*) x 

J — OC 

4(0)^(0), *t(*)F A (i)' 



Ho 



(3.21) 



Thereby both the charge and spin current across the con- 
tact can be calculated in terms of equilibrium correlation 
functions. 

To calculate these correlation functions, it is more con- 
venient to switch to a Euclidean Lagrangian approach. 
Note that we require correlators calculated not with re- 
spect to Ho, but to Hq. Thus we must consider 

Lfm = Lfm - dr dx ^FMpf , 

J J — oo 

Lqw = Lqw - J dr J dx LtiQW^V + h • 

Fortunately, this modification of the Lagrangians has no 
effect on the F and *Sf correlators. Physically, this is 
because the lead and wire correlators are calculated in 
equilibrium, so that the added terms act as potentials. 
They thus simply rigidly shift the spectrum of states on 
both sides of the contact, and those states raised or low- 
ered below the (equilibrium) chemical potential are filled 
or emptied, respectively. In general, this would induce 
some weak change in the correlators due to energy de- 
pendence of the density of states. For our model, how- 
ever, the correlators of interest are strictly unaffected. 
Formally, this follows since the transformations 



f(x) -> exp [ir z hfmx/v] f(x) , 
ip(x) -> exp it z ( hq W + h ■ ^ J ^ 



transform Lfm —> Lpu and Lqw — + Lqw, leave F 
and \& invariant, and respect the boundary conditions 
at x = 0. Note that calculating expectation values using 
these transformed fields (governed by Lfm and Lqw) in 
a functional integral formalism naturally produces corre- 
lators normal ordered with respect to the shifted fields. 
This captures correctly the physics of filling/emptying 
the shifted energy eigenstates discussed above. 

From the above discussion, it is apparent that the real- 
time correlator appearing in Eqs. (3.20) and (3.21) can 
be calculated using the pure, unpolarized Lagrangians, 



Lfm and Lqw corresponding to Eqs. (2.5) and (2.1), re- 
spectively. Their SU (2) invariance therefore implies 

( [4(0)^(0), (t)F A (t)l ) 0(-t) = 6 aX S^iC(-t) , 

\ L J / H o 

(3.22) 

where 

iC(t) = e(t)([B(t),B^(0)]) (3.23) 
is the standard retarded Green's function of the operator 
B = F T f * T . (3.24) 



The choice of spin components in Eq. (3.24) is arbitrary. 

Substituting Eq. ( j3.22| ) into Eqs. fl3.20| ) and fl3.2l| ), and 
using W(-t) = U(tj~we find 

/CO 
dtTr (U(t)W^Wa) C(t) , 
-oo 

/oo 
dtTr (U(t)W^W) C(t) . 
-oo 

It is helpful to express the matrices appearing in these 
expressions in terms of projection operators, 



s=±l 

u{t) = ex p 



, 2 u 

S | U'S 1 



s'=±l 



with v s = (1 + sh ■ a)/2 defined analogously to u s , see 
Eq. (2.6). Then it becomes straightforward to compute 
the averages 

Tr (u s v s >) =^(l + ss'rh ■ h^j , 

Tr (u s v s kt) = — (^sm + sh + iss'rh x h^j , 

and hence the tunneling spin current is 
G 



Ju 



■E 



(Prh + sh)Jm C(V + hs/2 + iS) 

(3.25) 



-Psrh xhRe C(V + hs/2 + iS) 

Similarly, the charge current is 

(I) = -G^2(l + Psm-h) lmC(V + hs/2 + iS) 



The quantities G and P were defined in Eq. ( |2.9| ), and 
we use the Fourier convention 



C(w) = / dt city 



(3.26) 
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The terms involving Im C are not surprising, since this 
is directly proportional to the spectral function of B, 
and hence has a simple interpretation in terms of tun- 
neling via Fermi's golden rule. For these terms, we can 
use well-known results discussed below. However, the 
terms involving the real part of C correspond to ex- 
change processes generated by tunneling, see Refs. U| 
and |22j. That tunneling indeed causes effective exchange 
couplings (even in the absence of a "bare" exchange 
coupling) follows already from the simple RG equations 
( |3.2[ ). As the physical effe cts of exchange are included via 
the boundary condition ( |3.16 ) with the 50(3) rotation 
matrix IZ^), we drop the terms proportional to ReC in 
the spin current ( 3.25 ). 

After some algebra, we then obtain the tunneling spin 
current as 



(./tun) — 



J^(Pm + sh)X a (V + hs/2,T), (3.27) 



For the case of a low-conductance conta ct, the exchange 
angle is small, $ <C 1, and Eq. ( |3.3C ) can be further 
simplified to 



J = — h x rh — ~Y^[ p m + sh]2 a (V + sh/2,T) . 

S 

(3.32) 



From Eq. (3.32), #/47r can recognized as acting as a sort 
of dimensionless spin conductance - proportional in fact 
to the " spin mixing con ductance" of Ref. The re- 
lations ( 3.31 ) and ( |3.32 ), together with the results in 
Sec. IV, provide the basis for our subsequent discussion 
of the setup in Fig. [|. 



IV. HYDRODYNAMIC DESCRIPTION OF BULK 
PROPERTIES 



where we have defined the function 



J a (U,T) = GlmC(U + iS) 



(3.28) 



Gk B T(k B T/D) a sinh Ir^l 



\2k B TJ 



r i 



Q 

2 



eU 



2nk B T I 



with the bandwidth D and the LL end-tunneling expo- 
nent a > 0. We note that in SWNTs D k, 1 eV, while 
a 1.1 according to the experiments reported in Ref. ||. 
In the limit eU <C k B T, the function (3.28) becomes 



T 



CJT 

yfeWr 2 (i + a /2) , 



and in the opposite limit, 
GU 



(eU/2nD) a 



The LL magnetization away from the contact can now 
be related to the spin chemical potential using the LL 
spin susceptibility x = 1 / (27r-u) , 



M=J R + J L = X h 



(3.29) 



Thereby we obtain the spin current J injected into the 
LL from the FM at any given contact for arbitrary ex- 
change coupling Omitting the expectation values for 
brevity, we find 



J=^-Sh+(l-S)J tnn , 



(3.30) 



where S = (1Z— 1 ) / (TZ+ 1 ) is a real antisymmetric matrix. 
Similarly, the injected charge current is 



\l + sPfn-h)l a (V + hs/2 1 T) . (3.31) 



In the previous section, we discussed how to determine 
the charge and spin currents in the neighborhood of a 
contact in terms of the local charge and spin chemical po- 
tentials fi and h of a LL. To complete the formulation of 
the full transport problem, we need to understand how to 
relate these quantities at different points within the LL. 
This is the subject of this section. In doing so, we assume 
that the length of the system is always long compared to 
some characteristic dephasing length beyond which the 
behavior is incoherent and hence classical. For a LL, we 
expect the dephasing length is set simply by the thermal 
length scale v/k B T. At low temperatures, this length is 
very long, and we thus are interested in the constitutive 
laws governing the LL on long length scales. 

On such long length scales, we expect a rather clas- 
sical description to apply both to charge and spin. For 
charge, classical behavior follows due to dephasing. For 
spin, more care must be taken, due to the non-commuting 
nature of the spin operators. On long length scales, how- 
ever, the total spin hs within any region is large. For 
large s 3> 1, all three components of the spin can in fact 
be specified with very good accuracy, the uncertainty 
being of 0(1/ s). Thus the long- wavelength "hydrody- 
namic" equations will involve simultaneously all three 
components of the magnetization and the corresponding 
currents. 



A. Operator equations of motion 

To construct the hydrodynamic equations, consider 
first the operator equations of motion for the spin cur- 
rents. These are obtained from the usual Heisenberg 
equations, dfJ^/iife) = i[H, Jyn(x)\, with the Sugawara 
form of the spin Hamiltonian.0 Writing 



magn i 
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the spin part of the LL Hamiltonian Hq is 



Ho= 2 I dx : Jr ' Jr + Jl ' Jl 



(4.1) 



The backscattering contribution has already been spec- 
ified in Eq. (2.2), and we now also include an external 
magnetic field B acting on the QW. For that purpose, as- 
suming a static and homogeneous field, we add the term 



Hmagn = ~fJ>B J dxB ■ M 



(4.2) 



with the LL magnetization (3.29). Here we have ab- 
sorbed the electron <?-factor g e into a renormalized Bohr 
magneton [Ib = g e e?i/2rnc. The equations of motion are 

(d t - vd x )J L = bvJ R xJ L - bvd x J R - hbJl* B , (4.3) 
(d t + vd x )J R = bv,J L x J R + bvd x J L - pb B J R x B . (4.4) 

Here terms like Jl x Jl are absent by virtue of normal- 
ordering. Their absence follows also by comparing the 
equa tions of motions obtained from the Sugawara form 
( |4.l| ) to the strictly equivalent equations of motion of the 
corresiX)nding SU(2) level k = 1 Wess-Zumino-Witten 
actionfl 



Taking the sum of Eqs. (4.3) and (4.4), we find 



where 



d t M + d x J = -hbM x B , 



J = (1 + b)v{J R - J L ) 



(4.5) 



(4.6) 



Comparing Eq. (4.6) to Eq. ( 2.12 ), we see that backscat- 
tering leads to a renormalization of the spin current. This 
is a "backflow" effect, s imi lar to those in Fermi liquid the- 
ory. For B = 0, Eq. (4.5) represents the standard spin 
continuity equation. Of course, the magnetic field will 
then spoil spin current conservation. 



Taking the difference of Eqs. (4.3) and (4.4) gives 



d t J + (1 - b 2 )v 2 d x M = bvM x J - fi B J x B . (4.7) 

For B = 0, we have a conserved spin current J but bulk 
precession of the magnetization around the fixed spin cur- 
rent. 



B. Hydrodynamics 



Note that Eqs. (4.5) and (4.7) are operator identities, 
not equations for the expectation values of these quan- 
tities. We will call these expectation values the classi- 
cal values. As argued at the beginning of this section, 
to describe the bulk physics on long length scales, the 
hydrodynamics should be phrased in terms of equations 
of motion for these classical variables. I n th e absence 
of backscattering, b = 0, Eqs. (O) and (4.7) are both 



Hbs = 



linear, and so taking their quantum expectation imme- 
diately gives the correct hydrodynamic description for 
the classical values. At zero temperature in the linear- 
response limit, moreover, the RG analysis demonstrating 
that b is marginally irrelevant implies that this b = 
hydrodynamics remains qualitatively correct. In general, 
even in this case there will also be finite renormalizations 
of physical quantities, but for b <C 1, these are expected 
to be small. 

At T > 0, and possibly also in non- linear response at 
T = 0, however, the hydrodynamic equations are gener- 
ally corrected by dissipative terms. Formally, these exist 
due to the fact that (M x J) £ (M) x (J) . Physically, dis- 
sipative corrections to Eq. ([O]) describe processes caused 
by the backscattering in which a non-zero spin current 
can decay. That the backscattering interaction mediates 
such processes can be seen by expressing it in terms of 
Fermions, 

—j J dx [Wr^rA^li + (TH.) 

+ A^rA^l]- (4-8) 

Consider an initial state containing one right-moving 
Fermion with up spin and one left-moving Fcrmion with 
down spin, thereby carrying a net spin current J z > 
but no magnetization. Acting on this state, i?b s flips the 
spins of right- and left-moving Fermions, thereby con- 
serving the magnetization but reversing the spin current 
J z — > — J 2 . Through a combination of such processes, it 
is natural to expect a finite lifetime for the decay of an 
initial spin current. 

Further, on physical grounds, we can express the rate 
of such decay processes based on Fermi's golden rule. In 
particular, we expect the inelastic spin current relaxation 
rate, I/t™, to be proportional to an expectation value 
that is quadratic in H^ s . More formally, we can deter- 
mine the lifetime using the leading non-constant correc- 
tion to the Fermion self-energy, which is the usual two- 
loop bubble in Fig. ^. Either way, since I/t™ is quadratic 
in 6, scaling determines the form of the decay rate, 



1/ T } n = Ab 2 k B T/h, 



(4.9) 



where we have used ksT to provide the energy scale 
needed from scaling. In addition, we have neglected 
more subtle logarithmic corrections expected on general 
groundsEJ The order unity numerical prefactor A is not 
obtained reliably by this simple argument, but a crude es- 
timate may be obtain ed u sing Fermi's golden rule. From 
the first term in Eq. ( |4.8| ), we may consider the rate for 
a single right-moving electron with down spin and mo- 
mentum k to flip it's spin, simultaneously creating a left- 
moving electron-hole pair: 



dkdqdq' , . 

— . . „ — {Air)6{vk — vk — vq + v q ) 
(27r) J 



x {2n)S{k + k' - q - q')f(vk') 
x [l-f(vq)}[l-f(vq')} , 
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where /(e) = \ j(e t l kBT + 1) is the Fermi function. This 
can be evaluated in th e lo w- energy limit vk <C k B T to 
give the result in Eq. (4.9), with however a surprisingly 
small prefactor A = 1/32tt. 

We proceed at th is st age on phenomenological grounds 
by modifying Eq. (4/7) by hand to include spin current 
relaxation, 

d t J+ (1 - b 2 )v 2 d x M = -J/tj + bvM x J- fijj/x B , 

(4.10) 

which provides a detailed (but approximate) descrip- 
tion of the crossover between ballistic and diffusive spin 
transport at low temperatures. In general, elastic impu- 
rity scattering processes will also relax the spin current, 
where in fact no spin-orbit interaction is required. This 
can be included according to Matthiessen's rule, 



1 



1 



1 



(4.11) 



Since the same impurity scattering processes relax both 
spin and charge currents, we expect 1/r} 1 to be compa- 
rable to the elastic scattering rate deduced from charge 
transport. 



Taken together, Eqs. (4.5) and (4.10) provide a start- 
ing point for the investigation of the spin hydrodynam- 
ics of a LL in the presence of electron-electron backscat- 
tering. A useful check is that these equations correctly 
recover the Landau-Lifshitz dynamicsEj of classical fer- 
romagnets at non-zero temperatures. Considering for 
simplicity B — 0, in the low-fre quen cy, linear response 
limit, the non-linearity in Eq. (4.10) is small, and the 
time-derivative of J can also be neglected. Then one can 
solve for J perturbatively in the non-linearity to obtain 
to leading order, 



a s = xD s = VTj/2n 



hv 



2nAb 2 k B T 



(4.15) 



Note that the spin conductivity has dimensions of length 
as expected in ID, and is essentially given by the mean- 
free-path for decay of spin currents, £ , = vt 7 . 



V. APPLICATIONS 

As an application of our general formalism, we now 
consider transport for a LL connected to two FM reser- 
voirs with magnetization directions rh\ and rh 2 for ap- 
plied voltage V, see Fig. [j]. The FM magnetization unit 
vectors mi :2 are tilted by the angle < 9 < it, so that 
rh\ ■ rh 2 = cos 9. For simplicity, we assume identical low- 
conductance contacts on both sides such that /iQw = 
and the exchange angle is small, t) < 1, Furthermore, 
we assume P < 1 for algebraic simplicity, although our 
results actually hold somewhat more generally.E3 The lat- 
ter condition is probably fulfilled in any practical appli- 
cation and ensures that the spin chemical potential is 
small, h/V -C 1. The validity of the latter condition is 
then self-consistently checked below. For h/V <C 1, we 
can use the following expansion for Eq. ( |3.28| ), 

l a {(V - sh)/2,T) = I a (V,T) - shG a (V,T) . (5.1) 

Here the current (per spin channel) for non-magnetic 
leads is 



I a (V,T)=I a (V/2,T) 



(5.2) 



such that for parallel FM magnetizations, 9 = 0, the 
charge current 1(0) = 2I a results in the absence of 
backscattering. Furthermore, the respective conductance 



J = -D s d x M - bvTjD s M x d x M , 
where the spin diffusion constant is 



hv 2 



Ab 2 k B T 



(4.12) 



(4.13) 



For SWNTs, £> <§; 1, and it is appropriate to approximate 
1 — b 2 w 1. In addition, the golden-rule estimate is at 
best valid for d « 1, and we therefore ignore the tiny 
renormalization v — > v' = (1 — b 2 )v. Inserting Eq. ( 4.12j ) 
into Eq. (4-5|) then indeed gives the usual Landau-Lifshitz 
equation.!^ 

We are predominantly interested in steady-state sit- 
uations in which both M a nd J are time- independent. 
Using Eqs. fl3.29|) and fl4.10D, we find 



d x h + —J= (b/v) (h + —^-—B ) J 
a s \ b 



(4.14) 



G a (V,T) = dI a /dV = (e 2 /2irh) 9o 



(5.3) 



Similarly we define a dimensionless contact conductance 
g via G = (e 2 /2-Kh)g. 

From Eq. ( |3~32; ), we can then write down the spin cur- 
rent Ji through the left contact, taken at x = with the 
local spin chemical potential h\ = h(0), and likewise J2 
through the right contact at x = L with I12 = h(L), 



1? - 

Ji = — h\ x rhi + PI a rhi — G a h\ 
4ir 



■h 



47T 



hi x m 2 + PI a m 2 + G a h 2 



(5.4) 
(5.5) 



The signs are chosen such that cur rents are oriented from 
left to right. Similarly, from Eq. ( 3.31 ), the charge cur- 
rent flowing through the device follows, 



where the linear-response spin conductivity is given by 
the Einstein relation 



1(9) 1 PG a% . 

= 1 ; All • Tfl\ 



2I a 



la 



(5.6) 
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where we have exploited current conservation, 
hi ■ rh\ + h% ■ rh% = . 



Next these relations describing the spin chemical po- 
tential and the spin current at the boundaries need to 
be related by virtue of the hydrodynamic description of 
Sec. [V. In the steady state, the basic hydrodynamic 



equations are Eq. (4.14) and 



d x J 



^hxB 
2nv 



(5.8) 



see Eq. (1.5). We first consider transport in zero mag- 
netic fie ld, a nd later on exte nd th e an alysis to finite fields 
in Sec. [V_b|. In Sees. [VA] and [Vj| the effects of the 
backscattering interaction are neglected, b — > 0, so that 
the spin resistivity vanishes. This is expected to be ap- 
propriate for short-to-intermediate length L. T he effects 
of b ^ will then be addressed in Sees. |V C| (without 
dissipation) and VD (including spin diffusion). 



A. Effects of Spin-Charge Separation 

We start with the simplest case of b = and B = 0, 



where the steady-state equations (4.14) and (5 



are 



solved by a constant magnetization and hence spin chem- 
ical potential, hi — h^ — h, and conserved spin current 



J. To compute the charge current (5.6), we then need 



to find h which in turn is determined from spin current 



conservation, Ji = J 2 . Using Eqs. (5.4) and (|5.5p, we 



then obtain three equations, namely Eq. (5.7) and 



4?r 



h ■ (mi x m,2) + G a h ■ (mi — 1112) 

-2PI a sm 2 (9/2) = , 
— cos 2 (9/2)h ■ (mi — m 2 ) — G a h ■ (mi x m.2) = . 

47T 

From these relations, with 1(0) = 2I a , the current results 
in the form 



m =1 _ p 2 tan 2 (9/2) 



1(0) ' ' tan 2 (9/2) + Y a 
Here the quantity Y a reads 

1? 



(5.9) 



W) = l + (- 



(5.10) 



For eV -C k B T <C D, using Eq. (5.3), this becomes 
(2d/gf 



Ya ~ r 4 (l + a/2) 
while for k B T < eV < D, 



(k B T/D)~ 



(1 + ay 



For a Fermi liquid, Yq = 1 + (2-d/g) 2 is related to the 
dimensionless spin mixing conductance 77 of Ref. f2G by 



(5-7) y = \rj\ 2 /Re(r]) , and Eq. (5.9) correctlsurecovers the 



current-voltage relation of a Fermi liquidEj in the limit 
a -> 0. 

In the interacting case, a > 0, however, for eV, k B T <C 
D, Eq. ( p.S| ) describes a drastically different behavior. 
Since Y a becomes very large for low energies, the spin ac- 
cumulation effect is completely destroyed for any 9 ^ tt. 
This remarkable phenomenon is entirely due to spin- 
charge separation, since only then the exchange coupling 
is so efficient at relaxing the injected polarized tunneling 
spin current. The absence of spin accumulation is then 
a direct signature of the presence of spin-charge sepa- 
ration and would allow to experimentally establish this 
phenomenon in a spin transport experiments Only for 
9 = 7r, one gets the standard 1 — P 2 suppression in the 
current. At low applied voltage, eV <C k B T, the jump in 
1 (9)/ '1(0) from unity for 9 < n down to 1 -P 2 at 9 = tt i s 



only smeared out by thermal fluctuations, see Eq. (5.10), 
and therefore becomes very sharp at low temperatures. 
Next we self-consistently check on the magnitude of 



h/V. Multiplying J\ - J 2 = by h yields with Eqs. ( ]5_J ) 
and (5.5), 



h 



PI a ( tan 2 (9/2) 
2VG a \tan 2 (9/2) +Y Q 



1/2 



(5.11) 



Note that M = x h = h/(2irv), so that Eq. fl5.1l| ) 
also describes the spin accumulation in the QW. Now 
I a /VG a equals 1 at high temperatures, and 1/(1 + a) for 
k B T < eV. Therefore h/V < P/2, and for P < 1, the. 
assumed smallness of h/V is self-consistently verificd.ESl 
Finally, we explicitly write down the spin current, 



- PI a , . s f , ^ tan 2 (61/2) 



(5.12) 



Note that from Eq. (5.7), this implies that h and J are 
orthogonal. 



B. Magnetic field dependence 

Next we consider a different and probably more fea- 
sible experimental setup aimed at revealing spin-charge 
separation. Instead of changing the angle 9 between the 
FM magnetizations, a simpler setup could work with a 
fixed angle 9 but employ the additional magnetic field 
B = BB. We assume that the bulk FM magnetizations 
mi, 2 are not affected by this magnetic field, and consider 
6 = 0. Under the se co nditio ns, t he steady-state hydro- 
dynamics in Eqs. (4.14) and (5.8) describes a precession 
of both right- and left-moving spin currents around B. 
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Therefore, when moving from the left to the right con- 
tact, the magnetic-field dependent precession phase 



7 = /isBL/hv 



(5.13) 



accumulates, so that j/2tt is essentially the ratio of the 
Zeeman energy to the level spacing. Since a field of 1 
Tesla corresponds to 0.058 meV, a sizeable precession 
phase can easily be achieved for tube lengths in the mi- 
cron range and standard magnetic field strengths. 

we can relate /12 and J2 at 



From Eqs. (4.14) and (5 



the right contact to the respective quantities at the left 
contact. Some algebra leads to 

I12 = COS7/11 + (1 — cos 7) (B ■ hi) B — sin 7 J\ x B , 

(5.14) 

and similarly 

J-2 = COS7 Ji + (1 — COS7) (B ■ Ji) B — sin'-fhi x B . 

(5.15) 

To compute the current, we then have to compute hi -rhi 
as outlined in Appendix ^| This is in general cumber- 
some and requires a numerical analysis. We thus restrict 
our attention to the special case 9 = tt, where the FM 
magnetizations rhi — — m 2 are antiparallel. Note that 
for 9 < tt, the spin accumulation effect is quenched any- 
ways by spin-charge separation, and then the magnetic 
field should have virtually no effect for a > 0. The case 
9 = tt is therefore also of most interest. 

First we consider a fixed field strength B. Let us sup- 
pose that B has been adjusted such that the precession 
phase is 7 = tt, but one allows for an arbitrary field di- 
rection B. In that case, the equation in Appendix |A| can 
be solved analytically, and the current is 



I(P)/2I a = 1 - P 2 



1 + Y a tan 2 ' 



(5.16) 



where Y a (V, T) is defined in Eq. ( 5. 10| ) and B-fhi = cos/3. 
If the magnetic field is parallel to the FM magnetizations, 
(3 — ► 0, a maximum spin accumulation effect is recovered, 
I(j3)/2I a = 1 — P 2 . Upon tilting the magnetic field, how- 
ever, since the quantity Y a diverges at low energies, spin 
accumulation is quenched again. This is very similar in 
effect to varying the angle 9 between the FM magneti- 
zations, but may be easier to implement experimentally. 
For (3 = 0, there is no precession since hi is parallel to 
mi, 2 and hence to B. Therefore spin accumulation is not 
affected by the magnetic field. From this observation, we 
then expect that the suppression of spin accumulation at 
(3 holds in fact for a broad regime of magnetic field 
strengths, and no special fine-tuning to 7 = tt should be 
necessary. The behavior of I((3) is qualitatively different 
depending on whether spin-charge separation is realized 
or not. 



Another possibility consists of fixing the field direction 
B and then measuring the I —V characteristics for differ- 
ent magnetic fiel d stre ngth B or, equivalently, precession 
phase 7, see Eq. ( 5.13 ). For simplicity, let us assume that 
B is adjusted perpendicular to rhi. Under this condition, 
the current can be found in closed form again, with the 
result 



2I„ 



1 



1 + F( 7 ) 



(5.17) 



where we use 



F( 7 ) 



(2V. 9a ) 2 -l + y a cos- 2 (7/2) 
( ffct /2^) 2 r Q +cot 2 ( 7 /2) 



(5.18) 



Note that Eq. (5.17) predicts a sharp negative transverse 
magnetoresistance peak near B = (and again at peri- 
odic intervals). In particular, at low temperatures, and in 
the linear bias regime, Y a , g~ x ^> 1, and one can estimate 
the full- width half-maximum (FWHM) of this feature in 
1(B) by 



(A 7 ) 



FWHM 



TT 



l) 

-in 



-1/2 



(5.19) 



which therefore decreases as T a as T — > 0. The height 
of the peak is, however, temperature independent. More 
generally, the current ( |5.17| ) is a periodic function of the 
precession phase 7, with f (7) = for 7 equal to an 
integer multiple of 2tt. In that case, the spin accumu- 
lation effect is maximal, and thus the current minimal, 
I( , y)/2I a = 1 — P 2 . On the other hand, for 7 being an 
odd multiple of tt, Ffr) — > 00 and thus a complete sup- 
pression of spin accumulation obtains. 

Remarkably, the way that one interpolates between 
these two limiting cases by varying 7 strongly depends on 
the correlation strength a, see Fig. ||. The 7 dependence 
of the current-voltage relation can therefore again reveal 
spin-charge separation. While for a — 0, the current is a 
smooth periodic function in 7, for a > and low ener- 
gies, the periodicity is turned into a series of very sharp 
dips at 7 equal to integer multiples of 2tt. This can be 
seen from Eq. (5.18), yielding for T — and eV D, 



F( 7 ) ~ V- 



-4a 



1 + (i9/4tt) 2 



! (7/2) 



(ii/47r) 2 +cot 2 (7/2) 



Therefore the spin accumulation effect is quenched un- 
less cot(7/2) diverges. Measuring the magnetic field dif- 
ference corresponding to the distance A7 = 2tt between 
two dips can also provide a di rect estimate of the con- 
tacted length L from Eq. ( 5. IS ) . 
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C. Backscattering: Dissipationless Precession 

For the remainder of this section, we address the con- 
sequences of a finite backscattering coupling 6, taking for 
simplicity B = 0. We start with the zero-temperature 
limit where the spin resistivity vanishes and Eq. (4.14) 
predicts bulk precession of h(x) around the fixed spin 
current, 



d x h = (b/v) h x J 



(5.20) 



Clearly, h ■ J is also conserved. It is then useful to intro- 
duce the quantities 



C b = 
Ctf — 



bgPLD 
2d 



(1 + a)g 



(5.21) 
(5.22) 



which provide dimensionless measures of the backscat- 
tering strength (C&) and of exchange (G&). The problem 
is then fully specified in terms of four dimensionless pa- 
rameters, namely the LL exponent a, the tilting angle 9, 
and of course Cb and C& . 

Similar to Sec. VB , precession of h(x) around the con- 
served spin current implies a precession phase when going 
from the left to the right contact, 



Aip = bJL/v 



(5.23) 



The spin chemical potentials hi and /12 can now be de- 
termined from spin conservation, J\ = J2, and the pre- 
cession equation (5.20). In addition, a self-consistency 
condition arises since J appears itself in the precession 
phase ( |5.23| ). This self-consistency implies that we are 
dealing with a nonlinear transport problem for b > and 
has far-reaching consequences. For convenience, use the 
abbreviation 



X 



hG 

Pin 



(5.24) 



with < X < 1, which provides a dimensionless mea- 
sure of the absolute value of the spin chemical potential. 
Then the charge current (5.6) is 



I = 2I Q (1 - P 2 X 2 



(5.25) 



From Eq. (5.4), the absolute value of the spin current can 
be expressed in terms of X as well, 



J 2 = (PI a ) 



($X 



(l-X 2 ) 



(5.26) 



where g a is given in Eq. ( |5.3| ) . In Appendix [B], we derive 
the following self-consistency equation determining the 
current, 



sin 2 (6>/2) = 



X 2 + {dX/2g a f 
1 + {$X/2g a y 



cos 2 (A<^/2) . (5.27) 



Using the precession phase (5.23) with Eq. (5.26), one 
can t hen s olve for X and thereby obtain the I — V rela- 
tion fl5.25p. 



Remarkably, solution of Eq. (5.27) predicts a multi- 
valued I — V relation. This is a direct consequence 
of the non-linearity of this spin transport problem in 
the presence of backscattering. Under an exact calcula- 
tion, we would in fact expect a unique answer for the 
curre nt, s ince thermal or quantal fluctuations around 



Eq. ( 5.27 ) should stabilize just one solution out of the 
multiple branches. Such a calculation could be performed 
in principle following a Keldysh approach or employing 
Langevin-type equations, but seems difficult to pursue in 
practice. A simpler approach based on an (approximate) 
free energy principle could obtain from the analogy to 
SNS junctions in Sec. VII A or the related "tilted wash- 
board" picture of Josephson junctions above the critical 
current. In the latter system, a free-energy principle is 
well-known to apply and to correctly resolve multistabjA-, 
ity questions posed by the equations of motion alone.EZI 
The difficulty with this approach for the spin transport 
problem at hand is to describe the dissipative terms in 
such a free energy. While the bulk energy is known, 
the boundary terms are less straightforward to handle, 
and, unfortunately, a rigorous resolution of this question 
has to remain open. Below we shall argue on intuitive 
grounds as to which of the multiple branches is realized. 

In ge neral , one needs to numerically find the solutions 
to Eq. ( [5.271 ), e.g. using a Newton- Raphson root-finding 
algorithm .E3 Numerical results accurately confirm an an- 
alytical solution possible for low energies, eV <C D, and 
a > 0. In the following, we focus o n th is regime of most 
interest, where Y a ^> 1, see Eq. ( [5.1C ), and search for 
solutions that fulfill also Y a X 2 ^> 1. Note that otherwise 
the effect of backscattering is negligible in any case. Un- 
der these conditions, the self-consistency equation (5.27) 



is solved by the precession phase taking only one of the 
following discrete values, 



Aip = (2n+ 1)tt - 9 



(5.28) 



where n — 0, 1, 2, ... is a winding number counting the 
number of full precession cycles of the steady-state bulk 
magnetization as one proceeds from the left to the right 
contact. In principle, there is also a set of so lutio ns ob- 
tained from the substitution 9 — > — 9 in Eq. ( 5.28 ) . For 
9 = and 9 = ir, both s ets co incide, and for < 9 < n 
we expect that only Eq. ( [5.28 ) gives stable solutions. 

Then the following general picture emerges. Focussing 
for concre tenes s on the case 9 = 0, the self-consistency 
equation (5.27) is solved either by X — or by A(p = 
(2n + 1)tt. As the voltage V is increased, first we have an 
arbitrary precession phase Aip < tt that increases with 
V. At the same time, Eq. (5.27) enforces X = 0, leading 
to the standard b — current I = 2I a . As the preces- 
sion phase hits Aip = n at the voltage V = Vq(9) (see 
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below) , the spin current J is locked at a fixed value such 
that Atp remains constant when further increasing the 
voltage. To keep J constant, however, the charge cur- 
rent I (or , equ ivalently, the quantity X) has to adjust 
from Eq. (5.26). Since this leads to a quadratic equation 
for /, there are two possible solutions for /. However, 
one of them would lead to unphysical currents exceeding 
21 ai and is disregarded in what follows. As the voltage is 
now increased up to V\ (9) , the precession phase Ay? = 3n 
becomes possible, and the above picture is re-iterated. 

For arbitrary 6, Eq. ( |5.2S| ) then predicts the following 
current- voltage relation. For V < Vo(9), the b = cur- 
rent (^1]) is realized. Upon increasing the voltage above 
the threshold Vq(9), however, sawtooth-like oscillations 
appear. In the window V n (9) < V < V n +i(9), we obtain 



D. Finite-temperature dynamics 



I(V) P 2 \ 1 (eV/4irD) 2a 

2I a 2 C$ 



(5.29) 



+ 



1 



(eV/4TrD) 2a 

1 /(2n + l)7T-0 
G& \ C b (eV/4nD) 



1/2- 



In the low-energy limit, the voltages V n (9) are given by 
eV n {9) (2n + l)7r-0 



4tt£> 



(5.30) 



and the I—V relation for V n (9) < V < V n +i(9) simplifies 
to 

T(y\ P 2 / , \ 

±1 = 1 - — (l + VI - (V n (0)/V) 2 ) . (5.31) 

Note that the V n (9) ~ 1/L and hence can in principle 
be made arbitrarily small simply by increasing the QW 
length L. Therefore the effects of backscattering become 
very important in sufficiently long QWs. Estimating the 
period AV — V n+ \ — V n corresponding to a full preces- 



sion cycle from Eq. ( 5.30 ) for typical SWNT parameters, 
we find AV w 10 to 100 mV. Measuring the oscillation 
period AV could provide useful information about the 
backscattering interactions and the exchange angle. 

The non-sinusoidal oscillatory I—V relation ( 5.31 ) 
is depicted in Fig. |^ for typical SWNT parameters and 
9 = 0, 7r. Apparently, backscattering has a dramatic in- 
fluence on spin transport not anticipated from thermo- 
dynamical considerations. To experimentally observe the 
predicted sawtooth-like oscillatory I — V relation, how- 
ever, it will probably be necessary to measure at very 
low temperatures using rather long and clean QWs. Note 
that these oscillatory behaviors are a non-equilibrium ef- 
fect not present in the temperature dependence of the 
linear conductance. 



As outlined in Sec. IV B| , for finite temperature, 
backscattering causes spin diffusion, and we now address 
the effects of spin diffusion on the current-voltage rela- 
tion, focussing on zero magnetic field, B = 0, where the 
spin is still conserved. The self-consistency equation in 
this case is derived in detail in Appendix ^|. To solve 
this lengthy equation, we again restrict ourselves to the 
low-energy regime, fc^T, eV <C D, with a > 0. 

For T <C T*, we then basically recover the results of 
while for T>T*, with the crossover tempera- 



Sec. |VC 
ture 



h v 



2ir 2 k B Ab 2 M 



1/L 



(5.32) 



a spin-diffusion-dominated regime emerges. This tem- 
perature is defined by n&L = a s (T*), see Appendix |^, 
with the spin conductivity ( 4.1 5| ) . In the following, we 
focus on the regime T T* , for which the I—V relation 
can be written as 



ijy) 

2/„ 



— (l + Vl - A n 



where 



2hG a f [(2n + l)ir - 9}ev 
b~PT n 



(5.33) 



(5.34) 



Remarkably, since a s ~ 1/& 2 , any dependence on the 
backscattering coupling b drops out in this temperature 
regime. In principle, the current is then again multi- 
valued and indexed by a winding number n. Here the 
appropriate threshold voltages, above which the respec- 
tive current can be realized, follow from the condition 
A n < 1. On physical grounds, in this spin-diffusion lim- 
ited transport regime, we expect that only the lowest 
winding number n = is realized, leading for V > Vo(9) 
to the current, 



^ = 1 - ^ (l + yfi (WF) , (5-35) 
where we assume kgT <C eV and use the voltage scale 



eVo(9) \ 



2+a 



(l + a)he 2 v 2 (l-0/Tr) 
2a s Lb 2 P 2 GD 2 



(5.36) 



Note that the right-hand-side of Eq. ( 5.3(: ) is propor- 
tional to (ksT/ D)(hv/ LD), with a prefactor of order 
unity. Hence this spin-diffusion limited regime should be 
accessible to experiments. 



VI. EXTENSIONS 

The approach developed in this paper is very flexible, 
and can straightforwardly be extended to describe a va- 
riety of other physically relevant situations. We sketch 
some of these extensions in this section. 
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A. Bulk contacts 

In many cases, particular in experiments on carbon 
nanotubes, contacts are made not to the ends of the quan- 
tum wire but to points in the "bulk" . Much of the pre- 
ceding theory applies to this case as well, but there are 
some additional complications. First, let us reconsider 
the problem of a single contact at x = x\, this time in 
the bulk (0 < x\ < L) rather than at the boundary. Be- 
cause there is no boundary condition relating right- and 
left-moving fields, there are more independent couplings. 
In general, the contact Hamiltonian H c has three dis- 
tinct contributions, neglecting redundant forward scat- 
tering terms which give spin-independent phase shifts, 



H r — H, 



H c 



H 



bsl 



Tun neling is described by the tunneling Hamiltonian 
( |2.7D , where \& = ip(xi), and we allow for different hop- 
ping matrix elements for right- and left-moving states 
of each spin polarization, tg~ . There are also two 
distinct exchange terms, 



a=R/L 



(6.1) 



Finally, there are local single-particle backscattering pro- 
cesses described by 



H hs i = Y Zs^rUs ■ aV L + h.. 



(6.2) 



where the projection operator u s is given in Eq. (2.6). 
These terms arise since the presence of a contact in- 
evitably leads to local disorder within the QW. 

For a non-interacting QW, all three contributions are 
on an equal footing, as they all involve Fermion bilin- 
ears. When interactions are present in the QW, how- 
ever, this is no longer the case. In fact, in the bound- 
ary RG framework, they scale completely differently. In 
particular, the tunneling Hamiltonian H tU n is irrelevant, 
the boundary exchange H cx is marginal, and the single- 
particle backscattering iibsi is relevant. Thus provided 
all terms in H c are comparable, i?bsi will dominate at 
low energies, where the distance to another contact, the 
inverse temperature, and any inverse voltage are all suf- 
ficiently large. The effect of such relevant backscattering 
terms has been studied extensivelyO We expect that the 
final result in the low energy limit is to completely sever 
the LL into two halves at x = xi. In that case, one 
can effectively ignore at very low energies half of the LL, 
namely the one not connected to a closed circuit, and 
treat the other half using the end-contact phenomenol- 
ogy of the previous sections. 

The latter discussion presumes that the £ s are substan- 
tial. In many cases, however, it is natural to expect that 
in fact |£ s | <C \K a \, |ii a ^|. For nanotubes, if the character- 
istic scale of the contacts is substantially larger than the 



inter-atomic dimensions of the nanotube, the £ s , involv- 
ing matrix elements that oscillate on the atomic scale, are 
considerably suppressed. Similarly, in contacts to semi- 
conductor quantum wires with widths large compared to 
the Fermi wavelength, £ s can be suppressed due to the 
smoothness of the effective potential at the contact. It 
is therefore of interest to describe the problem in the 
absence of single-particle backscattering, £ s = 0. The 



equation-of-motion methods of Sees. Ill and IV can then 



be straightforwardly extended to describe bulk contacts. 
For illustration purposes, consider in particular the case 
of multiple tunneling contacts at points < Xi < L, 
neglecting for simplicity bulk backscattering. Then the 
equations of motion for the chiral currents are 

(d t + avd x )J a = -fJ,BJ a x B (6.3) 

+ ^ 8(X - Xl ) \jT n {Xi) + Kafhi x jj , 



where a = L/R = qp, a nd th e tunneling current J^ un (xi) 
is determined from Eq. ( |3.27 ), with chemical potentials V 
and h approp riate to that contact and chirality. Further- 
more, in Eq. (3.27), one has to replace the end-tunneling 
exponent a by the bulk-tunneling exponent abuik < o^.Q 
Similarly, the chiral charge currents Ir/l obey the equa- 
tions of motion, 



(dt + av c d x )I a = y^6(a 



(6.4) 



where v c is the charge velocity. Each pair of equations 
can be combined to give equations for the spin and charge 
densities and currents, which can be solved in the steady 
state knowing the tunneling spin and charge currents at 
each contact (determined by the voltages) and the bound- 
ary conditions Ir = Ii and Jr = Jl at x — and 
x = L. Precession, diffusion, and spin-orbit scattering 
(see below) can also be included simply by adding the 
appr opri ate terms to the right-hand-sides of Eqs. (6.3) 
and (O). 



B. Nanotubes and Flavor 

While the general methodology and physical results 
applied above pertain to any interacting QW, some dif- 
ferences exist that should be taken into account in ap- 
plying the formalism in detail to nanotubes. In particu- 
lar, a "pristine" SWNT has not one but two ID bands 
crossing the Fermi energy, arising from the sublattice re- 
flection symmetry of the graphene lattice. Thus in fact 
the low-energy description of SWNTs requires an ad- 
ditional "flavor" index A = 1,2 on all electron fields, 
ipaa — > ipaAa- Moreover, the low-energy Hamiltonian de- 
scribing the nanotube in the absence of backscattering 
and magnetic fields respects the full chiral U(4) x U(4) 
symmetry of arbitrary separate U(4) rotations of right- 
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and left-movers in the combined spin-flavor space. This 
symmetry implies that, to a good approximation (since 
backscattering terms are weak), not only the SU(2) spin 
currents Jr/l but the full SU(4) spin-flavor currents. 



J, 



AaB0 _ . ,/,t 



(6.5) 



are conserved. In the ballistic limit, these currents satisfy 
chiral wave equations away from the contacts, 



(6.6) 



and a full solution of the transport problem in the steady 
state requires imposing constant values of each of these 
(16 + 16 = 32) chiral currents between contacts and/or 
ends of the nanotube. Moreover, SU(4) generalizations 
of the contact exchange H ex can be expected. 

As we have seen in the simpler single-channel case 
above, backscattering terms, even when weak, can lead 
to significant effects in long tubes. Because the backscat- 
tering interactions do not respect the (accidental) SU (4) 
symmetry but only the physical symmetries, there are in 
fact a variety of independent backscattering couplings. 
Thus, in general, the extension of i?bs to include flavor 
is rather complicated - the interested reader will find the 
(8 or 11, depending upon whether the nanotube is un- 
doped or doped) diverse backscattering terms enumer- 
ated in Refs. [Io| and ||o|. These interactions lead to a 
variety of generalized "precession" terms in the ope rato r 
equation of motion, which must be added to Eq. (3.6). 
These have the general form 



{d t Tvd x )Ji 



where the J^bcd ~ bv are related to the detailed form 
of the backscattering terms, and repeated indices are 
summed. 

In the hydrodynamic equations for the classical val- 
ues of these fields , we expect damping terms similar to 
that in Eq. ( 4.1 0| ) . In fact, there are two distinct sorts 
of damping processes which need be considered. First, 
like for the flavorless problem, backscattering terms lead 
to decay of the (non-chiral) SU(4) currents, J AaB P = 



AaB/3 



f 



(a) 



J. 



BDEF°R 



Aa) jCfj.B/3 jEfiFa 
JCAEF J R J L 



AaDu jEQFv 



v{ JA aB _ ?A«BP 



J? aDp ), which can be included via a lifetime 



tj. Second, however, unlike in the SU{2) case, the flavor 
densities themselves can decay via backscattering. Only 
the true spin magnetization, 



M = (4 
and charge density, 



'AaA/3 , fAaAfls -. /9 



n = -e{JT + J£ aAa ) , 

are requ ired to obey continuity equations, namely 
Eq. (4.5) and the usual charge continuity equation, by 



spin-rotational and U(l) invariance. The remaining or- 
thogonal linear combinations of the spin-flavor densi- 
ties can themselves decay, with some "flavor decay rate" 
i/ T / ~ 1/tj. For sufficiently long nanotubes and low 
energies, Vrf <C 1, these flavor densities become neg- 
ligible between contacts, and we expect to be able to 
simply ignore the flavor currents. Given the smallness of 
the backscattering couplings in nanotubes, however, this 
may occur only at very low voltages and temperatures, 
and for very long tubes. 

A proper treatment of these effects at intermediate 
length scales is technically rather complicated, and be- 
yond the scope of this paper. Nevertheless, the extension 
is in principle straightforwardly based on the techniques 
developed here. It is amusing to note that the issue of 
flavor currents is actually a concern even when all the 
leads are ordinary paramagnetic metals, since even such 
contacts generically do not respect flavor. 



C. Spin- flip scattering 

Up to this point, we have assumed that the QW in con- 
sideration is itself spin-rotationally invariant. In general, 
quantum mechanical spin-orbit coupling mixes spin and 
orbital angular momentum, leading to a small violation 
of spin conservation. This can occur both as a bulk and 
a boundary effect. For semiconductor QWs, the spin- 
orbit effects arc well understood. For SWNTs, they are 
expected to be extremely small. Indeed, in an ideal flat 
sheet of graphene, spin-orbit effects are negligible due 
to the high symmetry of the zone-boundary wavevector, 
the p z nature of the electronically active carbon orbitals, 
and the small atomic number of carbon. From a sim- 
ple tight-binding treatment (not shown here), we indeed 
find a vanishing effect for ideal graphene. In SWNTs, 
bulk spin-orbit effects can then in principle solely occur 
due to the curvature of the nanotube, to phonon distor- 
tions, and to defects in the tube. The first two factors 
are probably negligible as they are suppressed both by 
the smallness of the relativistic nature of spin-orbit cou- 
pling and by the smallness of the nanotube curvature and 
the electron-phonon coupling, respectively. Tube defects 
generally destroy the local symmetry of the lattice and 
allow some spin-orbit scattering. However, these same 
defects also elastically scatter ordinary momentum, so 
that one may estimate the spin-orbit scattering rate as 
1/t"so = e/^j 1 , where £«1 reflects the relativistic nature 
of the microscopic spin-orbit coupling. Since the elastic 
mean-free-path of SWNTs is known to be of the order of 
microns, the corresponding spin-orbit scattering length 
must be orders of magnitude larger, and hence also of 
no importance for current tubes which are at best a few 
tens of microns long. Parenthetically, we note that in 
multi-wall nanotubes, disorder is more important, and 
hence spin-orbit scattering may be of more relevance. In 
any case, it is straightforward to include the effect of spin 
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non-conservation theoretically by modifying Eq. (|4.5|) to 



8 t M + dj= -M/t so - n B M x B 



(6.7) 



Probably more significant is spin-flip scattering at the 
boundaries of the QW and contacts, which are often 
much more disordered than the bulk of the QW. Such 
processes can be incorporated by a renormalization of 
the contact parameters P, G, and z?, see Ref. Efl. 



D. Dipolar fields 

A unif orm magnetic field has already been included in 
Eq. (4.2). In general, magnetic contacts give rise to a 
spatially-varying dipolar field acting on the electron spin 
in the QW. Provided the variation of this field is smooth 
on the scale of the Fermi wavelength of the QW, however, 
the hydrodynamic treatment of the magnetic field can be 
applied to this case as well, simply letting B — > B(x) in 
Eqs. (4.5) and (4.10). Because the characteristic spatial 
scale of variation of the magnetic field is the size of the 
leads themselves, this condition should be amply satisfied 
for nanotubes and most semiconductor QWs. 

One can get an idea of the magnitude of the effect of 
dipolar fields by considering an idealized uniformly po- 
larized spherical FM of radius r c and magnetization rho 
end-contacting the QW at x = 0. The external mag- 
netic field of such a sphere is a pure magnetic dipole, 
hence B ~ |mo|/(l + x/r c ) 3 , with the usual dipolar de- 
pendence on the orientation of mo- To maximize the 
effects of the dipolar interactions, assume an Fe contact 
(which has larger magnetization than Co or Ni), and a 
radius r c comparable to the length L of the QW, so that 
B ~ \fho\ over the whole length. For low-temperature 
iron, \fno\ ~ 0.17 T, and thus Eq. ( 5.13| ) leads to the 
typical phase change, estimated for a SWNT, of 



Tlfpolar _ ^ 003 (L\ 



2tt 



(6.8) 



Thus the effect of dipolar fields is probably negligible. 



VII. DISCUSSION 

We conclude this paper by establishing a connec- 
tion to Andreev currents in superconductor-normal- 
superconductor (SNS) junctions and by pointing out 
some open questions. 



A. Analogies to Superconductors 

An interesting view on many of the above results fol- 
lows through an analogy to Andreev processes in ballistic 
SNS junctions. Consider for example the device depicted 



in Fig. [T[ Without loss of generality, we choose the mag- 
netization of the left FM lead as fci\ = x, and that of the 
right lead in the x — y plane, rhz — cos(#)i — sin(0)y. Ne- 
glecting electron tunneling, magnetic fields and backscat- 
tering, the Hamiltonian fully decouples into spin and 
charge components. Because the spin Hamiltonian is in- 
dependent of electron-electron interactions, we are free 
to model it using effectively non- interacting Fermions tp. 
Note that this in no way implies that the QW is non- 
interacting, but simply represents the physics of spin- 
charge separation. With boundary exchange couplings 
K\ and K 2 , we thus have 



-iv / dx 
Jo 
K2 

' 2 
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h.c. 



h.c. 



(7.1) 



Equation ([m]) must be supplemented by the boundary 
conditions ipR = ipL at x = and x = L. Consider now 
the (spin-down) particle-hole transformation, 



(7.2) 



which retains canonical anticommutators for ij) and pre- 
serves the boundary conditions. Under this transfor- 
mation, the kinetic terms in H a are invariant, but the 
boundary terms become anomalous, 



H a = -iv / dx^r z d x tp (7.3) 
Jo 

- (ViA(0) + K 2 e l0 K(L) + h.c.) , 

where the pair field is 
1 



(7.4) 



and we used the boundary conditions to remo ve t he fac- 
tor of 2 in the magnetic exchange. Equation (|7.3| ) is the 
Bogoliubov-deGennes Hamiltonian for an SNS junction 
in the limit of large normal reflection. The latter limit 
is implied by the boundary conditions V'i? = at the 
ends. 

The presence of the pair-field terms leads to Andreev 
reflection at the boundaries. As is well-known, such an 
SNS junction carries an equilibrium current I for any_£ 
which is not a multiple of 2ir. In particular, we expectEil 
for 9 of order one, 



I = ev ( il^T z %j) 



K 1 K 2 } ev 



(7.5) 



in the fully coherent limit, eV, ksT <C v/L. To translate 
this result back into the spin problem, we can make a 
dictionary relating quantities in the two pictures. Some 
interesting variables in the spin problem are 
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J z = ^W z cr z i> = I/2e , 
m z — \^o z i\j = - : := n/2 , 



We see that J z and m z correspond to the charge cur- 
rent and density, respectively, in the transformed vari- 
ables. Thus the FM-LL-FM device indeed carries a non- 
vanishing z-axis spin current. Note, however, that the 
in-plane magnetization corresponds to strange "large- 
momentum" pair fields in the analog SNS system. 

For comparison to the results of the previous sec- 
tions, note that our hydrodynamic treatment gives zero 
spin current at zero ap plied voltage. This is not in- 
consistent, because Eq. ( |7.5| ), which is exact in equilib- 
rium, predicts a spin current J z ~ v/L that vanishes 
in the thermodynamic limit. More precisely, the hy- 
drodynamic results require incoherent transport, which 
holds e.g. for eV 3> v/L. The hydrodynamic approach 
predicts J z ~ GPV, whi ch c an be crudely matched to 
the "Andreev" prediction ( |7.5| ). In particular, the hydro- 
dynamic and "Andreev" currents are comparable when 
eV ~ v/PL, that is essentially at the boundary between 
the coherent and incoherent regimes. One learns from the 
SNS mapping that the spin current is actually enhanced 
by coherence. 

In an SNS junction, one expects a proximity-effect in- 
duced pair field within the normal region. Naively one 
might therefore expect some uniform bulk magnetization 
in the x — y plane. This conclusion is, however, false, as 
can be seen by rewriting the superconducting pair field, 



(7.6) 



The pair field thus maps back to the 2kp oscillatory com- 
ponent of the x — y magnetization, A 
but not to the uniform one. 



L 2k f 



B. Outlook and open questions 

Let us finally summarize some of the open questions 
from our point of view, and provide an outlook. 

One rather obvious concern might be the incoherent 
nature of transport assumed in our study. For sufficiently 
long quantum wires and/or low conductance of the con- 
tacts, it certainly is appropriate to assume a two-step 
sequential transport mechanism through the FM-LL-FM 
device. What happens if one has coherent transport? 
The latter situation could arise for higher-transparency 
contacts or at very low energies. However, from the anal- 
ogy to SNS junctions, we expect that our main conclu- 
sions are qualitatively unaffected by coherence, and we 
therefore do not expect a dramatic change. Nevertheless, 
it would be interesting to study this question in detail. 



For non-interacting electrons, this could be done in the 
framework of a Landauer-type approach, see also Refs. |2(] 
and [2^. 

A related issue concerns the role of charging effectsEl 
where transport through the LL is hindered by Coulomb 
blockade. For low-transparency contacts, these effects 
are known to be crucial at energy scales below the charg- 
ing energy E c , which can be estimated for a SWNT of 
length L, radius R, and background dielectric constant 



E c = (e 2 /nL)\n(L/R) 



and typically is of the order of a few meV. Charging ef- 
fects are washed out by intermediate-to-high tempera- 
tures, and could in principle be avoided altogether by 
using higher-transparency contacts and/or long tubes. 
We note that ^larging effects also tend to destroy spin 
accumulation,^ and therefore one has to be careful that 
they are not present when experimentally testing for 
spin-charge separation. However, since they manifest 
themselves through quite pronounced dependencies on 
external gate voltages, this issue is not expected to cre- 
ate serious difficulties in practice. Furthermore, although 
our present theory does not include charging effects, this 
could be accounted for easily via a proper treatment of 
the zem-modes in the bosonized version of the Luttingcr 
liquidilH 

A very interesting extension of the methods of this 
paper is to problems involving mesoscopic ferromag- 
netic contacts which are sufficiently small so that their 
magnetization becomes dynamical. Here the quantum 
wire/nanotube would mediate an effective RKKY-type 
interaction between the FM magnetizations, and novel 
transport phenomena can be anticipated. 

It would now clearly be of great interest to experimen- 
tally study the scenario put forward here. The probably 
best candidates for such experiments are single-wall car- 
bon nanotubes, which should offer the unique possibility 
of observing spin-charge separation directly on a single 
ID quantum wire. In addition, the effects of backscatter- 
ing were shown to imply rather dramatic consequences 
for spin transport, such as a sawtooth-like oscillatory 
current-voltage relation. Such spectacular consequences 
of the electron-electron interactions have not been pre- 
dicted previously, but should be observable for long nan- 
otubes at very low temperatures. 

Future work should also address in detail the 2D gen- 
eralization of these ideas, which seems particularly inter- 
esting in the context of some theories of high-T c super- 
conductivity. We hope that our paper has convinced the 
reader that spin transport in strongly correlated meso- 
scopic systems represents an exciting area of research 
that leads to both fundamental insights and technologi- 
cally useful devices. 
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APPENDIX A: TRANSPORT IN A MAGNETIC 
FIELD 

Here we outline the main step in the derivation of 
the I — V characteristics in a magnetic field for 6 = 0, 



Using the abbreviations (5.24), W = J/PI a > 0, and 
Z = $/ (4irG a ), the relat ion h ■ J\ = gives h ■ rhi = X. 
Furthermore, from Eqs. (5.4) and fl5.5p, 



mi • Ji = to 2 • J 2 = (1 -X 2 )/W . (B3) 
Next we use that h ■ J 2 = implies 

h ■ rh,2 = XZ sin( Aip)m ■ J — X cos(Aip) , 
where m ■ J is determined by Eq. (B3). In addition, 



?7J2 • Ji = rh ■ J, see Eq. (B3), gives 



see Sec. VB. To do so, we first eliminate J\ 2 and h 2 



from spin current conservation using the spin currents in 



Eqs. (5.4) and (5.5), and the relations (5.14) and (5.15) 



We are then left with the following bulky relation deter- 
mining hi, 

§ — $ r, » - 

— COS7 (mi + 7712) x h± + — (1 — C0S7) (B ■ hijirhi x B) 

47T 47T L 

-{hi ■ (mi x B)}B] + (1 + G 2 a )smjhi x B 

+ (i?/47r) 2 sin7 (B ■ hi)(rhi x m 2 ) + (B ■ rhi)(m 2 x hi) 

i9 r- r - - 

+ — G„ sin7 hi{B ■ mi + m 2 )} — (B ■ hArhi 
4tt L 



-(hi ■ rh 2 )B 



2G n 



COS7/1! + (1 — cos7)(-B • hi)B 



h ■ (mi x m 2 ) = (2/XZ)[sm 2 {e/2) - X 2 cos 2 (A^/2)] 
+ X sin(A(/j)m • J . 

Finally, we may employ the relation rhi-J 2 = rh-J, which 
yields 

X Z cos(Aip)h ■ (mi x rh 2 ) 
— X Z sm(A(p)(m ■ J) h ■ (fhi — rh 2 ) 
-2sin 2 (6>/2) + 2X 2 cos 2 (A<^/2) (B4) 
— Xsin(A(^)mi ■ (h X J) = . 

Alternatively, one could use spin conservation of J - (mi x 
777,2), which produces the same answer. When simplifying 
Eq. (IB 4ft, it is helpful to use the following relation, 



= P/ Q |cos7mi — 7772 + (1 — cos7)(£? • rhi)B 

7? - -1 

H sin 7 [(_B • 7712)7771 — cos 6* + G a sin 7 mi x B > . 

47T J 



mi • (h x J) = - X 2 )/T4^ = -Zlrii • J , 

which follows from Ji • J = J. Straightforward alge bra 
then leads to the self-consistency equation Eq. ( 5.27 ) . 



This equation is then analyzed for 



7r and special 



choices for B in Sec. VE 



APPENDIX C: SPIN DIFFUSION 



APPENDIX B: DISSIPATIONLESS PRECESSION 

The algebraic manipul ation s necess ary t o obtain the 
self-consistency equation (5.27) in Sec. VC are provided 
in this appendix. With {hi — h 2 ) ■ J = 0, we can write 



hi-h 2 = l-(h- J) 2 cos(A(^) + (h- J) 



(Bl) 



Since hr Ji = h 2 -J 2 is conserved, by multiplying Eq. (5.4) 



by hi and Eq. (5.5) by h 2 , and exploiting Eq. (5.7), we 
get hi t 2 ■ J — 0. This in turn implies directly that we can 
relate h 2 to hi. With h = hi, we obtain 



cos(Aip)h — sm(Acp)h x j 



(B2) 



The unknown spin chemical potential h can then be ob- 
tained from spin current c onse rvatio n, J i = J 2 , with the 
currents specified in Eqs. (5.4) and (5.5). 



In this appendix, the technical steps in the derivation 
of the I — V chara cteristics in the presence of spin diffu- 
sion, see Sec. |VD|, a re given. The charge current can be 
written as Eq. (|5.25| ), but with a modified quantity X , 



v2 G a J* 

A = — ; hi ■ mi 



Pla 



(CI) 



In additio n, we use Ct — PI a L/a s , with the spin con- 
ductivity (4.15). To compute the current, we first need 
to express h 2 in terms of hi = h via the steady state 
diffusion-precession relation (4.14). For symmetry rea- 
sons, h 2 2 — hi 2 , since we consider identical contacts. 
That directly implies from Eqs. (|]J) and (5.5) that 



h- J =-h 2 -j = WC T /2 



(C2) 



where we use again W = J/PI a . With the precession 
phase Aip defined in Eq. ( |5.23 ) we then get instead of 
Eq. (H), 
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h 2 = cos(A(p)h - aia(Aip)h x J — WC T cos 2 (Aip/2)J . 

(C3) 

Combined with Eq. (|C^), this allows to express h 2 in 
terms of X 2 alone, 



(hG a /PI a ) 2 = X 2 - W 2 G a C T /2PI a 
Here J\ 2 = J 2 yields with Z = $/irG a , 
(l-X 2 )(l + Z 2 X 2 ) 



w 2 = 



l + G a C T (l + Z 2 )/(2PI a ) 



(C4) 



(C5) 



Then Eq. flBj) still holds. 

We now employ spin current conservation to obtain 
a closed nonlinear self-consistency equation for finding 
X 2 and thereby the current-voltage relation. With Q 2 = 
(G a /PI a )h-rh,2, the relation mi • J 2 = PI a (l — X 2 ) gives 
after some massaging 

X 2 - ZW( lf T sin(A^) - 4 sin 2 (9/2) sin 2 (A^/2) 

cos(A</?) + — (1 - X 2 ) sin(Av3) 



G u Ct 
PI a 
G u Ct 

PIa 



cos 2 (Av?/2)(l 



cos 2 (Av?/2) cos(0)(l + Z 2 X 2 ) = 



The second relation allowing to eliminate Q 2 comes from 
h- J 2 = WCt J/2, and reads 



1 



ZWG a C T 
2PL 



cot(Ay>/2) 



Q 2 + cos(Aip)X 2 



--sin(A^)(l-X 2 )A 2 



G a Cx 

PI a 
ZWG a 

2PI a 
G a Cr 



G a CT 

2PL, 



cos 2 (A<p/2)W 2 



cot(Ay>/2) X 2 + {2sin 2 (A^/2) 



cos 2 (A^/2)}(l-X 2 ) 



. 



PIa 

Eliminating Q 2 from these two relations gives the self- 
consistency equation for X 2 for arbitary temperature 
and applied voltage. The solu tions to this equation give 
directly the current via Eq. ( 5.25 ). One checks easily 
t hat t his reproduces the T — self-consistency equation 

We shall now evaluate the self-consistency equation 
in the spin-diffusion dominated regime characterized by 
T > T*, with the scale T* defined in Eq. flsll^ . This 
temperature results from f(T) — ZGq,Ct/ PI a — 1 for 
T = T*. Since f(T) > 1 for T > T*, the above equa- 
tions can be drastically simplified in this regime. 

For T ^S> T* , the self-consis tency equation is again 
solved by the discrete values (5.2S|) of the precession 



phase Aip indexed by the win ding number n. Since un- 
der these conditions, from Eq. flC5| ), the precession phase 
can be written as 



it is then straightforward to derive Eq. (5.33) in Sec. VD 
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FIG. 1. Proposed experimental setup (schematic). An individual SWNT or quantum wire is connected via low-conductance 
contacts to two ferromagnetic reservoirs, and I — V curves should be measured either in an additional magnetic field or for 
various angles between the magnetization directions in the leads. 




FIG. 2. Self-energy diagram to estimate the decay time for the spin current. 
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V/V 

FIG. 4. Current-voltage relation in the presence of backscattering for T = and 6 = 0, 7T. The scales are set by To = 2/q, 
and V = V (6 = 0). 
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